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SUMMARY 


j 

The  dynamical  ec^uatlons  of  motion  axe  formulated  for  a  body  duriug 
re-entry  into  the  atmosphere.  With  the  assumption  that  the  translational 
■ind  rotational  degrees  of  freedom  may  be  treated  independently,  an  analysis 
is  undertaken  of  the  oscillatory  motion  of  a  re-entry  vehicle  vhich  is 
spinning  about  its  longitudinEil  axis.  The  influence  of  the  spin  rate  end 
the  static  margin  upon  the  stability  of  the  transient  solution  of  the 
eqxiations  of  motion  is  considered,  and  it  is  concluded  that  satisfactory 
dynamical  behavior  of  the  re-entry  body  may  be  anticipated  if/ 
the  static  margin  of  the  vehicle  is  less  than  zero^ 

^  the  spin  rate  is  in  the  range  0.5  to  2  radians  per  second^  t 
the  lift  curve  has  a  slope  equal  to  or  greater  than  zero. 

An  examination  of  the  precessional  and  nutatlonal  motion  of  the  veMcle  is 
also  made,  and  the  effect  of  initial  conditions  upon  these  modes  of 
oscillation  is  indicated,  influence  of  spin  rate  and  static  margin 

variations  upon  the  fore  utions  of  the  equations  of  motion  Is  con¬ 
sidered,  and  the  following  {conclusions  are  drawn: 


o  spin  rates  on  the  c|:dor  of  0.5  to  2  radians  per  second  effectively 
'average  out'  re-entry-body  asyimaetries . 
o  for  a  spin  rate  of  one  radian  per  second  the  lateral  displacement 
of  the  re-entry  body  from  the  original  trajectory  plane  is  con¬ 
siderably  less  than  a  thousand  feet  at  i:iq>act  if  the  static  margin 
is  of  the  order  of  minus  six  inches  or  loss. 

Finally,  a  brief  analysis  is  made  of  the  planar  re-entry  case,  and  the 
results  are  coo^pared  to  those  obtained  by  previous  investigators* 
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dynanlc&l  equations  of  ckotion  are  foxmilatod.  for  a  body  dorlJOf 
re-entry  into  the  atzoospbere.  With  the  assumption  that  the  trzinsXatlonal 
and  rotational  degrees  of  fmodaa  aay  be  treated  Independttitly,  an  analysis 
la  luidertaken  of  the  oscillatory  notion  of  a  re-entry  vahicle  which  Is 
spinning  about  Its  longitudinal  axis.  The  Influence  of  the  spin  rate  end 
the  static  aargin  upon  the  stability  of  the  transient  solution  of  the 
equations  of  motion  is  considered,  and  it  Is  concluded  that  satisfactory 
dynaaical  behavior  of  the  re-entry  body  aay  be  anticipated  if: 

0  the  static  margin  of  the  vehicle  is  less  than  zero, 
o  the  spin  rate  Is  in  the  range  0.3  to  2  radians  per  second, 
o  the  lift  curve  has  a  slope  equal  to  or  greater  than  zero. 

An  exsmin&tlcsa  of  the  precesalcxial  and  nutatloml  ootlon  of  the  vehicle  is 
also  made,  and  the  effect  of  initial  conditions  upon  these  modes  of 
oscillation  is  indicated.  The  Inflxjoace  of  spin  rate  and  static  margin 
variations  upon  the  forced  solutions  of  the  equations  of  notion  is  con¬ 
sidered,  and  the  following  conclusions  are  drawn: 

o  spin  rates  on  the  order  of  0.3  to  2  radians  per  second  effectively 
'average  out'  re-entry-body  asynsaetries . 
o  for  a  spin  rate  of  one  radian  per  seccoid  the  lateral  displacement 
of  the  re-entry  body  from  the  original  trajectory  plane  is  con¬ 
siderably  less  than  a  th<xmaQd  feet  at  lxq)act  if  the  static  margin 
is  of  the  order  of  oinvui  six  inches  or  lees, 
finally,  a  brief  analysis  is  made  of  the  planar  re-entry  case,  axid  the 
results  are  coopared  to  those  obtained  by  previous  investigators. 
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^  ••  angles  relating  the  position  of  the  velocity 

-  vector  vith  respect  to  the  inertial  reference  axes 

-  angle  betveen  tb^  reference  plane  and  the  horiiontal  plane 


0  -  the  steady-state  fllgiri~p«±h  angle 

8S 


an^pOLar  rates  vith  respect  to  an  liMrtial  reference;  l/sec 
aboot  the  re-entry  body  orthogonal  axes 


Coordinate  Ax&e 


X  - 


T>- 
Z  - 

’j 


reference  axes  fixed  vith  respect  to  the  earth  and  considered 
to  be  inertial 


re-entry  body  axes  of  syasastry 


X  - 


ro-«itry  bwly  central  principal  axes 


Miscellaneotts  Terras 


g  -  acceleration  of  gravity  (ft/ sec  ) 
h  -  altitude  (ft) 

k  -  atxaosidtiere  density  exponential  (a  negative  nuiaber) 


V  -  re-entry-body  speed  (ft/sec) 


R  -  cpln  rate  paraaeter  equal  to 

X 


I  O 

X  X 
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I.  UJTBOIXJCTIOH 


One  of  the  problems  vhlch  axlse*  during  the  re<-eatry  into  the  stuios- 
phere  of  e  body  ia  the  dispersion  of  the  trajectory  due  to  aerodysauaic  lift 
forces.  Such  umranted  forces  Bight  result  frcoi  asysaetries  in  nass  distri¬ 
bution,  asywetriee  of  the  body  about  the  longitudinal  axis,  or  other  caxises. 
Since  very  kbsH  lift  forces  can  cause  errors  at  ispeunt  on  the  order  of 
several  nllesf^^  It  ia  desirable  to  eliainate,  or  at  least  to  reduce,  the 
effect  of  these  forces  »q)on  the  re-entry  trajectory. 

case  B»ans  of  doing  this  is  to  reduce  the  dispersion  of  the  re-«itry 
vehicle  by  causing  the  lift  vector  to  preceas  about  the  velocity  vector. 

The  lift  and  drag  forces  exert  torqxies  upon  the  re-entry  vehicle,  which, 
if  the  body  is  spun  about  the  longitudinal  axis,  will  cause  precession 
about  the  velocity  vector.  A  high  precession  rate  tends  to  "average 
out*  the  lift  rector,  thus  reducing  dispersion. 

Although  precession  of  the  lift  vector  about  the  velocity  vector  soeas 
to  offer  a  solution  to  the  problem  of  dlsx>erslon  due  to  aarodynaadc  forces, 
the  dynsaBle  stability  of  the  re-entry  body  is  influenced  by  the  spin  rate 
and  thus  should  be  investigated.  Furtheraaore ,  it  may  be  anticipated  that 
spinning  the  re-entry  body  will  introdxice  a  unidirectional  side  force  which 
will  cause  a  lateral  dlspiacement  of  the  vehicle.  This  effect  is  slailax 
to  that  which  occurs  in  the  case  of  projectiles  fired  from  rifled  guns. 
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n.  csRiYAyios  of  m  zqgmosa  of  mohon 


Figure  1  defiofts  tii«  eo<»‘dlzuit«  syvteaa  asd  goM  of  tho  razdables 
pertiaant  to  tb«  analysis  of  tha  rotational  sotion  of  tha  iplnning  ro^ozztry 
body  Tbs  X  axis  in  tbo  referencs  plana  is  eillgned  in  tb«  direction  of  tbs 


viBloclty  vector  at  the  'beglimini;  of  the  prohlen  and  Is  considered  to  be  an 
Inertisd  reference.  The  orleirtation  angles,  V,  0  and  relate  the  position 
ocf  the  re-entry  body  orthogoial  axes,  x  ,  y  and  z  ,  to  the  inertial  reference 

AS  S 

axes,  Z,  Y  ax^  Z.  The  angles  ^  and  specify  the  position  of  the  velocity 
vector  vith  respect  to  the  inertial  reference  vhile  the  angles  ^  &dcL  3  relate 


the  position  of  the  velocity  vector  to  the  re-entry  body  axes. 


Fig  I—  Re-entry  geometry 
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A  casjplffte  analyslB  of  the  re-entry  probleni  voold  re<i:atre  the  solutloa 
of  six  coxgjled  equatioafl.  However,  If  it  ie  asstnoed  that  the  three 
rotational  degrees  of  freedoB  have  a  negiigible  effect  upon  the  translational 
BWtlon  of  the  center  of  gravity  of  the  vehicle,  then  two  sets  of  three 
equations  each  say  be  solved  independently.  Solotioiis  to  the  flight  path 
equations  must  be  obtained  before  solving  the  rotatioiuU.  equations  of 
action. 


'Die  three  translational  and  the  three  rotational  equations  of  action 
are,  respectively: 

bV  «  Bg  sin  (6  +  0  )  -  D 


(1) 


=  Bg  005  (0^  +  0^)  - 


(2) 


hV^  cos  0^  =  -  (3) 

\diere  D  is  the  drag,  L  ,  the  lift  due  to  L-  the  lift  due  to  P,  and  0  the 

p  o 

angle  between  the  reference  and  horizontal  plane. 


I  -I  M 

\  .  z  y  X 

^  ^  I  ^  r 

X  X 


(M 


I  -I 


M 

JL 

I 

y 


(5) 


I  -I 


M 


0J„ 


y  X  .  z 

+  -‘‘v —  CJ  (J  =  v~ 

I  y  I 

z  z 


(6) 


where  M  ,  M  and  M  are  the  torques  about  the  x,  y  and  z  axes  respectively. 
X  y  z 


The  relationship  between  body  angular  rates  and  the  rate  of  change  of 


the  orientation  angles  my  be  deduced  by  an  exaaination  of  Fig.  1; 


5 


-  t 


-  Y  sin  0 


(7) 


Q  cos  f  4-  Y  sin  f  cos  0 


(8) 


*  Y  cos  t  cos  Q  -  0  sin  (9) 

With  the  &ssu^ptioa  thst  the  aa^^s  of  attach  aisd  side  slip  ejre  8»all,  the 

folXotving  eagpressions  SAy  be  found  bj  an  inspection  of  Fig.  1: 

-  «  +  0  (10) 

P 

0  »  (^  +  t)  cos  (ll) 

Th©  form  of  Eq.  (4^-6)  indicates  that  the  re-entry  body  axes  are 
principal  axes.  In  Pig.  2,  the  relationahip  between  the  principal  axes  and 
the  eaces  of  8>iasaetry  of  the  vehicle  is  indicated.  With  a  noaunlfona  mass 


V  c  os  . 


S  . 


Fig  2--Re-en^ry  body  axes 

distribution,  the  origin  of  the  principal  exes,  00,  E»y  be  displaced  frosa 

the  longltudii^  axis,  x  ,  while  the  x,  t  axes  may  be  rotated  through  an 

s 

angle  £  with  respect  to  the  x  ,  z  ,  axes.  Moreover,  it  may  be  expected 

s  s 

that  the  mosaents  of  inertia  about  the  y  and  £  axes  vxnild  not  be  e<tual. 

However,  if  6  and  £  are  oaiall  in  miiignltude,  their  products  with  other  ^ 

£ 
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quantities,  such  as ^  and  p,  laay  be  neglected,  and  the  difference,  for  a 
sjEaetrical  body,  betveen  and  becoaes  negligible. 

The  torques  abotft  the  x,  y  and  s  axes  are  thus 


H 

X 

I 

X 


(12) 


M  A  5 

^  r  *  "x  r  -  S.S 


(13) 


A 
_3 

y  I 


=  =  F 


(IM 


where  A^  is  the  static  aargin  of  stability  and  6^  is  the  dlapiaceaant  of 
the  center  of  mass  of  the  re-entry  body  frocj  the  longitudinal  axis 
syaoetry.  The  paraaeter  is  related  to  the  stability  derivatives, 


and  C 


The  coBsponents  of  lift  and  drag  along  x,  y  and  z  are 


F 


F  L_  cos  4  -  sin  »  -  Doc  sin  t  +  D  P  cos  t 

y  0  ^ 

F  -  -  L-  sin  if  -  cos  v  -  Da.  cos  »  -  D  0  sin  t  DG 
z  p 

Upon  assuming  an  exponential  atmosphere,  ,  L  and  D  have  the  following 

a.  0 


form- 


(15) 

(It) 

(17) 


(16) 


kh 


-  e  V  0 


(19) 


kn  I 
D  ^  e  V 


(20) 
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v^ro  k  is  a  negatirs  nuaber. 

Tha  pamwtars  asd  Bay  be  vrltt«n 

\  °  2  %  ^ref  (^) 
aod 


(21) 


^  "  2  S)  ^ref . 


(22) 


vh«re  la  the  lift-curre  alope,  la  the  drag  coefficient,  is  the 

aea-Xerel  air  deneity  atid  A  .  la  a  reference  area. 

rcr 

The  pax^icular  aet  of  rariables  utilised  to  re-entry  body 

stability  la  at  the  diecretion  of  the  investigator.  Thus  one  sight  use  cJ^ 

end  G3  ,  0  «ad  T,  or  ^  and  By  coaeiderlng  Kq^.  (2)  and  (5),  in  conjunction 
s 

vlth  Eq.  (>*20),  expressions  a»y  be  obtained  for^  and 

CoMbinlag  Eq.  (2),  (10)  and  (l8)  and  differentiating  yields 


^  —  V  e^i  +(tat 

81 


La 


-Of 


B8 


Utilising  Eq.  (5),  (6),  (8)  and  (9),  ve  find  for  0 


..  M  M 

0  -  cos  t  -  jp  sia  ♦  - 


I-I 


CJ. 


(23) 


(2^) 


Equation  (2^)  any  be  slaplified  by  the  substitution  of  Eq.  (13)  and  (14), 
in  eonjunetion  vith  Eq.  (8)  and  (9)  and  Eq.  (1^-20).  Ihus 


o’  «  ~  (V^)^  "  hi 


PC,  +  ^  7  e^tfC  -  0 
c  '^ss 


r 


-  / 


n  - 


I-I 

cj 

I  X 


Z  X  1 


)  cos  t 


(25) 


8 


Hovev«r,  frcao  Bq.  (3)  and  (ll) 


t 


m 


(2c) 


■Ria  introduction  of  Eq.  (25)  and  (26)  into  Sq.  (25)  leada  to  a  aocond -order 
equation  in  oC,  vith  6  and  0  coupling  tems. 

In  a  Banner  identical  to  that  indicated  above,  the  corresponding 
equation  in  0  my  be  obtained.  Thus 


CL  +cC 


ri  kh 


—  V  e  *(jJ 
B  n 


(27) 


.  (,-E^) 


6  ♦  —  T  3 
m 


■  *  O  ® 

t-  0  *0  +<:j  (t — i~)(“f - )  cos  f 

M  '^88  n  ^3- 


. 

dA„kh  „\„kh  2 

e  +  0 

5"  ^  *  *  Si 

-^3 

-j—  (— —  V  e  /  ♦  Kw  —  Ye  •  CJ 
dt  'a  'mb  n 

(28) 


(t-R  ) 
z 


5 

Z  X 


L.  'Ml'  f-r— >  ’ 


vhere 


I-I 

o  X 

R  --  -= — cj 
z  I  X 


0  -  ,4  cos  0 

88  V  '^88 

In  deriving  Eq.  (7')  and  (26),  the  folioving  assuaaptloas  have  been  Bade 
(l)  the  angles  CL.  0,  ^p’  ^  and  / are  ' small  angles';  (2)  the  aerodynaxalc 
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paraartera  aaid  are  tine  -  Invar  l«int ;  tb*  EOcaeirt  of  inertia  aixnxt  y  ift 
o<2ual  to  the  acBont  of  inertia  about  z,  the  oscillatory  notion  of  the  vehicle 
has  X30  effect  upon  V  or  0  ,  the  steady-state  fli^t-path  angle. 

6  B 

It  nay  be  noted  that  for  the  case  of  a  aca-splnDlng  re-entry  body,  Eq 
(27)  reduces  to  the  planar  for®  which  has  been  treated  by  H  Julian  Allen 
and  others . ^ 

Equaticna  (27)  and  (26)  are  coupled,  and,  for  the  case  in  vdiich  C»^  la 

equad  to  a  constant,  linear  with  variable  coefficients.  The  spin  rate, 

is  consteurt  if  the  torque  about  the  x  axle  is  zero,  since,  with  equal  to 

I  ,  -UJ  cross-coupling  due  to  CJ  and  cJ  occurs.  Frca  Sq  (12)  and  (l.  )  it 
Z  J  z 

can  be  seen  that  a  diaplacecant  of  the  CG  froa  the  axis  of  synsaatry  gives 

rise  to  a  torque,  M  .  However,  aaeuialng  that  5  is  a  quantity,  the 

X  z 

tortiuee  which  arise  frosa  such  an  effect  nay  be  safely  neglected  s.nce  terms 

of  6  eC or  6  &  ars  second  order.  Purtheraore  M  is  an  oscillatory  function 
Z  Z  X 

which  on  the  average,  approaches  zero. 

By  coabining  Eq.  (V),  (10)  and  (11)  we  find  for 


d  -  S  -  p<?t*  e  0  -  B  0  B  0  0"  *  0  (:v) 

t  P  2  P  ■  (  (  P 

« 

If  ve  again  neglect  hlgh-order  tenas,  v  Is  equal  to 

In  Appendix  A,  the  hoaogeueous  solutions  of  Eq  ( }  and  (26  )  are 
considered,  while  the  forced  solutions  are  treated  in  Appendix  Appendix 
C  ccaslders  the  aarodyaaalc  tenas  and  their  relative  magnitudes  for  a  t}-p1;su 
re-entry  body. 

A  planar  re-entry  analysis  Is  examined  in  Appendix  D,  and  the  results 
are  ccaspared  with  those  of  i.ir«vlous  studies. 

Qiere  is  one  special  case  of  interest  when  it  is  necessary  to  retain 
at  Least  the  secorid -order  terms  of  Sq.  (29)  That  is  when  is  equal  tc 

toro  Under  this  conditicsn  EJq.  (2^'  and  (26'>  are  coupled  non-.^lnear  equi-  ions 
with  tine -dependent  coefficients. 


\ 


III.  DISOJSSIOH  OF  THE  REOLTS 
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The  Transient  Solutlcuu 

The  ho«ogeneouB  solution*  of  Eq..  (27)  *nd  (28)  as  indicated  in  Appendix 


A  have  the  following  fora; 
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If  th«  forerltxxSy  of  the  re-entry  vehicle  is  conical,  then  the  aerodynaaic 


paraaeters  of  £q..  (30)  and  (31)  Bay  be  written  as  follows. 


d  »  y  e 


kh 


2N 


2 


(cot  B-1)  +  3511121  g  I 


A  ^  a  p 
+  cot^  B 


(a) 


CJ 


^o  *  ®  ^2  ^  ^  kh 


(b) 


p^e  p 

(cot  B-l)  V  e 


kh 


(c)  (32) 


"l 


I  O 

X  X 


21 


-  ^o  ®  ®  ^2  „2  kh  ^x  H: 


k  I 


dt 


(d) 


I  CJ 

X  X 

—  Ti 


A^  a  g 


I 


2N 


4 


dt 


(e) 


As  it  is  noted  in  Appendix  A,  the  two  exponential  dashing  terms  are  very 


nearly  identical,  except  when  the  static  margin,  approaches  zero.  151X13 
for  most  cases  the  stability  of  the  transient  solution  isay  be  deduced  by  an 
exaainatlcQ  of  the  tern 


d  dt 


"o 

175 

Lii  4  ___ 

The  integral  of  the  daaping  term  xsay  be  approximated  as  follows: 


t 


d  dt 


„  2  A  ^  a  cot^  B 

(cot^  B-l)  +  S - - = - 

36TriN  sin^  B  ^ 


( 


kh 

e  dh 


sin 


ss 


(33) 


since 


h  ■  -  y  sin  0 


36 


See  Appendix  C 
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Dova  to  altitudes  of  ^,000  feet  or  so,  the  steady-state  flight-path  angle 
changes  very  slowly.  Thus 


t 

■^5 


d  dt 


4  H  k  sin  6  !_ 

*  A  a 


88 


(cot  B-l)  + 


2  2 
A  a  cot  B 

X 


56TriN  sin  B 


kh 


kh 

0\ 

e  ) 


(34) 


'Where  sin  p  Is  an  average  quantity  over  the  altitude  Interval  considered. 

88 

since  k  Is  negatl-ve,  -the  as^Utude  of  'the  exponential  tern  vUl  decrease  vlth 
decreasing  altitude  as  long  eis  the  quantity  vithln  'the  bracket  Is  positive. 
For  cone  half -angles  of  4^  deg  or  less,  this  condition  Is  alvi^  fulfilled. 


Plguzre  3  Is  a  plot  of  the  noraallaed  eagponentlal  ejqxsnent  as  a  function  of 
the  cone  half -angle  for  a  particular  re-entry  body  configuration.  An 
exaalnatlon  of  Fig.  3  Indicates  that  for  a  vehicle  vith  a  A  of  -  0.5  feet, 
zero  dacplng  occurs  for  a  cone  half -angle  of  55  <ieg* 


Fig.  3 


-  Normalized  damping  function  vs 
re-entry  cone  half-angle 


RM-1863-1 

14 


!Ibe  roHLlnlng  factor  affecting  the  ax^Iitude  of  the  oscillatory  action 


is  the  ter»  l/ 


R  2 

X 


n  lA  2 

'  .  Frc*  Xq..  (32b ),  it  can  be  seen  that  CJ 


TJ  -  - 

approaches  zero  at  very  high  altitudes.  Thus,  since  the  initial  conditicsis 

2 

are  estsbliahed  vhen  CJ  is  approxljaately  zero,  the  aarlgum  valxie  attained 


r  «  R  ^  1/4 

by  the  tern  1/  IcOj^  +  -^  is  Va/R, 


1/2 


As  a  result  of  this  feurb,  even 


vith  zero  dsspl.ng,  the  aaxljnw  amplitudes  of  ^  or  0  never  exceed  the  initial 

value  of  the  envelope  of  the  oscillatory  action.  Although  the  aagnitude 
2 

of  CO^  decreases  after  the  peak  in  dynamic  pressure  has  been  passed,  it 
never  becomes  zero  at  altitxides  belov  this  point.  Figure  4,  a  plot  of 
nonaalized  frequency  function  versus  altitude,  illustrates  this  point. 

For  the  special  case  in  ^dnieh  the  static  laargln  is  zero,  Rq.  (50 )  and 
(31)  simplify  to  the  following  forms: 


k  sin 
<<-  e 


. ,  kh 

') 


(b"*  -  e  0) 


kh 


ffC  SI  h  sin  0 

~  sin  R  t  e 
R^  X  o 


SB 


and 


p  n  -  e 


k  sin  i 
'^88 


VV,  ^ 

/  kh  o 
(e  -  e 


) 

^  n  k  sin  ^ 

^  cos  R  t  .  g-  . 

X  X 


(s'*  -  e  “) 


88 


(35) 


(36) 


2 

fv. 


There  are  now  two  distinct  exponexxtial  terms  which  influence  the  stability 
of  the  solutions.  Sie  exponent  containing  A^  arises  due  to  C  ,  and  thus 

q 

this  term  is  always  convergent.  The  second  exponential,  however,  is  a 
function  of  the  lift  curve  slope,  C.  ,  which,  for  a  conical  body  with  cone 

half -angles  greater  than  45  deg,  can  have  a  negative  value  at  hypersonic 


normolized  frequency  function 
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With  a  eoae  half -angle  greater  than  k3  deg,  Sq.  (35)  indicates  that 
^haa  a  divergent,  noa-osciUatory  solution.  Of  course,  if  is  positive, 
the  re-entry  body  is  subjected  to  a  de-stahilizing  sooent,  and  divergence 
vill  occur. 

Before  considering  the  initial  conditions  of  the  problem,  it  is 

desirable  to  select  a  sequence  of  events  in  the  operation  of  the  attitude 

control  systea  of  the  re-entry  body.  Let  \is  assume  that  after  the  last 

stage  of  propulsion  has  ceased,  the  vehicle  is  rotated  by  reaction  forces 

until  the  longitudinal  axis  is  aligned  vith  the  direction  the  velocity 

vector  will  have  i^on  re-entering  the  atmosphere.  After  the  vehicle  is 

stabilixed  in  this  position,  an  angular  acceleration  is  iaparted  until  a 

* 

predetermined  angular  rate  about  the  longitudinal  axis  is  achieved.  At 
this  point  the  control  period  is  terminated. 

The  initial  conditions  that  have  been  selected  in  the  study  of  the 
transient  behavior  of  OC  and  B  are : 

attvtjorhsh 
o'  o 

cCa  qC  I  P  B  0 
o 

a 

♦  ♦  • 

CLm  ^  I  B  =  0 

o 

With  an  ideal  attitude  control  system,  all  of  the  fnltial  conditions  vould 
be  zero,  and  only  the  forced  solutions  of  end  B  vould  exist.  For  an 
actual  c<nitrol  system,  rate  and  position  errors  in  both  OL  and  B  should  be 
esQ^cted.  However,  since  the  primary  motion  of  the  longitudinal  axis  of 
the  re-entry  body  during  the  control  period  is  in  the  plane  containing 

The  dynsaics  of  the  spin-up  have  xiot  been  Included  in  this  report. 
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^  and  3  should  he  small  in  magnitvide.  and  may  be  neglected, 
o  o 

If  the  residueLL  body  rate  Is  very  saall  cooipared  to  the  x>osition  error, 
then  Eq.  (30)  and  (Sl)  may  be  approximated  as  follows; 


uid 


vhere  it  has  been  assxaned  that  R^>  >0. 

The  total  transient  angle  of  attack  is  thus  approximately 


Equation  (39)  indicates  that  vith  an  initial  of  zero,  the  lift  vector 
precesses  about  the  velocity  vector,  but  nutation  does  not  occur.  The  rate 
of  precession  may  be  fo\md  by  differentiating  the  angle  yhoee  tangent  is 
OC  Thus 

P  =  — 5 

+  3^  L  J 


(40) 
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OtillzlJig  Eq.  (37),  (3S),  aad  (39),  ve  find  for  p 

.  .  R  I  I 

p »  =  -  -|  +y^  +  -f 

0 

?rcm  Eq.  (4l)  It  can  be  seen  that  as  gtoes  to  zero,  p  approaches  zero.  As 

• 

A^  it  increased  in  a  negative  sense,  p  approaches  for  a  fixed  value  of 

R  .  The  cotq^llng  paraaster  R  ,  ^diich  is  equal  to  I  is  prinarily  a 

functioi  of  the  spin  rate,^.  Previous  studies  have  indicated  that  a  spin 

rate  of  ftroB  0.5  to  1  rad/sec  is  re(piired  to  reduce  the  effect  of  steady- 

state  trl»  anglas  of  attack  irtiidi  arise  due  to  various  1;ypes  of  re-entry- 

body  asyasetries.  Thus'  a  range  of  values  for  of  frosa  0.5  to  2  rad/sec 

vould  seas  adequate  from  such  considerations.  With  A^  equal  to  a-tnus  O.5 

feet  or  larger,  and  with  R^  in  the  range  of  values  indicated  above,  is 

very  large  coapared  to  R^.  Under  such  conditions,  in  the  altitude  region 

in  \diich  the  peak  dynsoic  pressure  occurs, O  ,  and  thus  p,  might  be  on  the 

a 

order  of  30  rad/sec.  A  precession  rate  of  this  magnitude  is  more  than 

sufficient  to  average  out  the  transient  lift  vector. 

Figures  5,  6  ai^i  7  arc  plots  of  ^  ^  ®'  ^^®^ton  of  altitude  for 

varlovis  coabinations  of  static  xaargins  and  nose-cone  half -angles.  The 

effect  of  the  spin  rate,d)  ,  upon  A,  JoL  is  Indicated  in  Fig.  8. 

»  ^,p  0 

If  the  body  angular  rate  at  the  end  of  the  control  period  is  such  that 

0 

is  not  negligible,  the  lift  vector  vtU  nutate  as  veil  as  process  about 
the  velocity  vector.  In  order  to  slaqplify  the  problem  the  altitudes  considered 
vill  be  restricted  to  200,000  feet  and  belcw.  With  much  larger  than  R^, 

Sq.  (30)  and  (51)  become: 


50  too  150  200  250  300  350  400 

Altitude, h  (thousonds  of  ft) 

Fig  5  -Total  transient  angle  of  attack  as  a  function  of  altitude 

Static  margin  =-0.5  ft 


6l 

T-£9eT-KH 


50  100  150  200  250  300  350  400  450 

Altitude, h  (thousands  of  ft) 

Fig. 6  — Total  transient  angle  of  attack  as  a  function  of  altitude 

Static  margin=  —  I  ft 
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X 

■If 


2,V 

"n  *T 


1/4 


and 


d  dt 


1/2 


cos 


Z  - 


R  t 

X 


ys  a. 


172 


cos  Z  sin 


R  t 

—  (42) 


•?/ 


d  dt 


3 


2 


21 


1/4  i 


7^ 


1/2 


sin 


Z  - 


R  t 

X 


y^  a. 


R  t 

Sin  Z  sin  V 


where  h  ^  200,000  ft  and 


/  "■> 


dt 


7 


kb 

“S 


k  s  in  j5 


ss 


(43) 


The  total  angle  of  attack  is 


-  I  /  4  dt 

,1/2  ®  o 


a,e 


and  the  precession  rate  is 

-n-®* 
n  -5- 


1  + 


■*  a  ® 


0  .,.2  V 


4  Q-^  a^t  B^t 

g-g-  sin  COB  ^  tta 

*  «  *x  °o 


1/2 


«  20-.  4cj  /rr  \2  n  R  "fc 

(<n  -  \  )  -  — sin  Z  cos  Z  -  -?[-§(  ~  )  sin^ 

“  a  R2\Q/  2 

o  X  \  0/ 


(44) 


(.45) 


4a. 


R  t  R  t  /dA2  «  R  t 

X  .  ^  tin  .f.  CPS -I- .  Sln^.f. 

X  o  \  o/ 


An  examination  of  Ec.  (44)  and  (45)  indicates  that  nutation  becomes  important 


2k 


vhcu  fi  very  large  sagaltudej  or  ^dieii  is  taall «  !Dis  exact  na'ture 

of  the  precessional  ootloa  depejads  upon  the  relative  asagnitudes  of  <3^^  R^ 

finA  cj  ,  It  be  seen  fron  Eq.  (^5)  that  p  nay  actually  be  negative  over 
Q 

portions  of  the  period.  When  is  zero,  p  equals  xainus  Th\x«  the  time 
history  of  the  laotion  of  the  longitudinal  axis  about  the  velocity  vector 
is,  in  geiMral,  very  coi^lex. 

The  envelopes  of  the  and  aiiniBUfi  sa^tlitudits  of  the  oscillations 

of  A  -  are  plotted  as  functions  of  the  altitiide  in  ?ig.  9*  the  special 

a,p 

case  in  ^diich  a  is  finite,  d  is  zero,  and  R  is  equal  to  one,  the  envelope 
O  O  X 

of  1^/2^  “  *  ftawtlon  of  altitude  is  given  by  Fig.  5-7. 

la  utilizing  Fig.  ^9,  it  should  be  kept  in  adnd  that  the  small  angle 

appxxnciisations  that  were  in  the  derivation  of  the  equations  of  motion 

in  Section  II,  break  down  at  altitudes  on  the  order  of  30, (XX)  to  50,CXX)  ft. 

An  examination  of  Pig.  5-9  indicates  that  increasing  the  spin  rate 

increases  the  aamlitxjde  of  A  A  coDg>arison  of  Fig.  5  of  Section  II  and 

a,S 

Fig.  10  of  Appendix  D  illustrates  the  differences  between  the  spinning  and 
non-spinning  cases.  In  the  limiting  case  as  R^  approaches  an  Infinite  value 
ve  see  from  Eq.  (30)  and  (31)  that 


m  k  sin  ^ 


kh 

/  kh  o  \ 

(e  -  e  ) 


ss 


'a. 3 


Thus  for  veiy  large  values  of  R  (or  conversely  for  very  low  values  of  ca  ) 
the  stability  of  the  oscillations  is  dependent  only  upon  the  lift  curve 
slope 


a 


The  normal  acceleration  load  acting  on  the  re-entry  body  is 


P..  8 


a, 3  2N 


kh 

e 


{kt) 


(<H'i 
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Bren  for  a  very  large  value  of  the  Initial  angle  of  attack,  sxich  as  20  deg, 
the  EaxlauB  value  of  n  for  a  typical  re-entry  body  varies  betveea  rcro  and 
five.  Zero  occuro  for  B  equal  to  1^5  deg. 

The  Forced  Solutions 

From  Bq.  (B-7)  and  (B-8),  ve  find  for  the  steady-state  solutions 


•  •  « 

^SS  Si  ^8S  /  Sd  > 

n  n 


(6  ) 

— —j - cos  O  t 

A  X 

X 


8 


ss 


R  ^  A_  (6  +eA  ) 

-  -  ^ o  t 


CO 


X 


Equations  (kS)  and  (49)  are  good  approximations  at  altitudes  of  250,000  ft 
and  less,  and  if  is  greater  than  sero.  At  very  high  altitudes  where  the 
aerodyxiaaiic  tenas  of  Eq.  (27)  and  (28)  approach  zero,  or  when  A^  is  zero. 


the  forced  solutions  are : 


aC  »  e 
88 


dt  ^  dt 

ic  d  0 


88 


P  «  0 
88 

Let  us  first  consider  Eq.  (46)  and  (49).  The  tenas  containing  5  are 

z 

those  that  arise  due  to  asyiacetries  in  the  mass  distribution  of  the  re-entry 
vehicle.  In  both  equations  the  terms  are  multiplied  by  sinusoidal  functions 
which  depend  on  the  spin  rate,  CJ^.  Reference  1  indicates  that  with  6^ 
ranging  from  O.5  to  1  rad/ sec,  the  miss  distance  of  the  vehicle  at  impact 


is  approximately  one  per  cent  of  the  value  which  occurs  without  spin.  The 
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tvo  tenna  of  Kq.  (kS)  axe  dvie  to  th«  lag  betveen  the  engiilar 
motion  of  the  velocity  vector  and  the  longitudinal  axis  of  the  re-entjry 
body. 

The  firat  term  of  Eq.  (J^9)  arises  due  to  a  cotqpllng  betveen  the  spinning 
body  and  the  aerodynamic  and  gravltatiaial  forces  acting  on  the  body.  It 
is  a  unidirectional  effect  ^ich  vlth  a  non-zero  lift -curve  slope  vill 


cause  the  vehicle  to  depart  from  the  original  trajectory  plane. 
Revrltlng  Eq.  {k$) 

-kh 

I  CJ  g  cos  0  c 

V  at'  '  o  • 

3 


(5) 


88 


(6  +t^  ) 

^  Sin  O  t 


b2) 


From  Sq.  (52)  it  can  be  seen  that  the  sign  of  the  unidirectional  term  depends 
on  the  sense  of  the  spin,  and  of  the  static  mrgin,  La  the  case  of 
a  projectile  fired  from  a  gun  vith  right-hand  rifling,  the  direction  of 
departure  is  to  the  right  of  the  trajectory  plane,  since  is  positive  for 
a  bullet.  A  re-entry  body  vlth  a  positive  spin  vwild  nove  to  the 
left  if  it  is  statically  stable  and  has  a  positive  lift-curve  slope.  The 
approximate  magnitude  of  the  effect  may  be  fotmd  as  follovs: 


g’  ^  .  ^2  kh  ^xSc  ^ss  \ 

Integrating  Eq.  53  yields  the  rate  of  departure  from  the  original  trajectory 
plane .  Thus 

1  u.  \ 


(K) 


8  . 


Hie  lateral  displacement,  as  a  function  of  time  is: 


(^) 
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Initially^  is  16  deg,  ^rixiie  at  iaqpact  0^^  is  90  An  order  of 

magnitude  value  of  S  may  be  obtained  by  approxliaating  (0  -  0  )>  since 

66  0 

the  Integral  of  Eq,.  (55)  nsay  be  expressed  as  follows: 


t 

r  dt .  t 

For  a  typical  re-entry  body,  with  H  equal  to  140,  a  mass  of  100  slugs,  and 

2 

I  equal  to  160  slug-ft  ,  ve  find  for  S 

R  (1  -  tan^  B) 

S  =  JS - ^ -  (77.5) 

X 

With  a  60  deg  cone  half -angle,  an  of  2  and  a  of  minus  0.5  ft,  8  is  a 
positive  620  ft.  It  should  be  reaambered  that  Eq.  (57)  is  uot  valid  when 
approaches  zero.  However,  from  an  examination  of  Eq.  (B-5)  and  (B-b)  it 
is  evident  that  may  becotas  very  small  before  the  approximation  is  no 
longer  applicable.  Ihus  the  xmidirectional  sideslip  term  can  load  to  large 
values  of  lateral  dispersion. 

When  is  zero,  Eq.  (50)  aai  (51)  represent  exact  solutions  of  the 

steady-state  equations  of  motion.  The  only  forced  disturbance  with  a  static 

margin  of  zero  is  in  (X  . 

8  8 

Prom  a  coB^jartson  of  Eq.  (49)  and  (51),  it  appears  that  as  A  becaass 
®aall  the  magnitude  of  the  unidirectional  term  increases  \mtil  a  maxisium 
is  reached.  As  A  continues  to  decrease,  B  also  decreases  until  it  Is 

X  ss 

zero  with  A^  equal  to  zero. 
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IV.  COKCLUSIOHS 


Frca  the  analysis  ot  the  transient  and  forced  solutions  for  oC  and  0, 
it  is  possible  to  draw  certain  conclusions  concerning  the  dynanics  of  the 
re-entry  vehicle's  angular  action: 

1.  Unless  the  static  icargin  is  positive  the  envelope  of  the 
oscillations  of  a  re-entry  vehicle  whoso  forebody  is  conical 
will  never  diverge  if  the  cone  half -angle,  B,  is  U5  deg  or 
less.  For  large  negative  values  of  the  static  aargln,  B  may 

be  larger  than  deg  and  the  envelope  will  still  be  convergent 
(sea  Pig.  7).  Ifaless  the  parameter  R  is  very  large,  or  the 

X 

static  margin  is  very  a*nall,  spinning  the  vehicle  has  no  effect 
upon  the  dynamic  stability  of  the  body. 

2.  The  amplltvide  of  the  oscillations  increases  as  R  increases, 

X 

At  altitudes  on  the  order  of  200,000  ft  or  less,  the  frequency 
of  the  oscillations  is  completely  determined  by  the  magnitude 
of  the  static  moment  unless  is  vary  large  or  is  very  small. 

3.  The  effect  of  the  forced  or  steady-state  solutions  of  and  3 


due  to  mass  asyraoetries  is  reduced  by  increasing  (1^,  since  spin 
tends  to  'average  ovrt'  the  lift  vector.  However,  the 
unidirectional  aideelip  term  increases  as  R  increases. 

X 

A  value  of  R^  of  frcm  O.t'  to  2  rad/ sec  is  adequate  from  the  viewpoint, 
of  averaging  out  re-entry-body  asyEoetrles.  With  in  the  above  range, 
the  unidirectional  aideallp  term  is  of  little  irportance  unless  the 
static  margin,  is  very  small.  However,  tne  problem  of  dispersion  becomes 
relatively  unimportant  if  the  re-entry  vehicle  Is  designed  so  that  the 
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llii;~cnrv«  slope,  ,  is  zero. 

Spia  aboxit  the  longitadlnal  axis  is  desirsible,  even  vith  C.  equal  to 

^(C 

zero,  in  order  to  prevent  a  large  initial  value  of  cC  frcai  arising  due  to 


saall  residual  rate  errors  at  the  end  of  the  control  period. 


Ihider  all  oircuxstances  it  is  necessary  that  be  less  than  zero. 

Free  Kq,  (37}  and  (38)  it  can  he  seen  that  a  of  zero  does  not  iJB^jly 
divergent  csclUaticms .  However,  the  daaq>lng  of  the  Biotion  is  very  weak, 
and  as  a  conseq^sence  relatively  large  values  of  ^  or  8  Bight  occur  in  the 
IcFver  atBOsphere.  IMs  vtmld  be  undesirable  due  to  the  problems  associated 


vith  serodynaaio  beating. 
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APPENDIX  A 

The  HcaaogeneoujB  Solutions  of  the  Bgxiatlaas  of  Motion 


With  the  assuEption  that  the  rotational  laotion  of  the  re-entry  body 
does  not  affect  the  time  history  of  V,  h  and  0,  these  variables  laay  be 
detenained  from  the  three  flight-path  equations  of  motion.  Squatioos  (2?) 
and  (28)  are  thus  cotqpled  linear  equations,  the  coefficients  of  ^diich  are 
known  functions  of  time. 

Revrritlng  Eq.  (27)  and  (28)  in  the  homogeneous  form 


where 


CC  ^  ^0^  +j^(X+  R 


3  +  d  p  +j^.  3  -  R 


[3  +  e  ^ 
1^4-  c^ 


(A-1) 

(A-2) 


_2  ..r  d  A  kh,  ^  ■'S.  „  lih 

R  =  -  R 

X  I  z 

\  „  kh 
c  a  —  V  e 
m 

As  an  example  of  the  method  oi  solution  that  is  to  bo  easployed,  the 
planar  oscillation  of  the  re-entry  body  witho'xt  spin  is  of  interest.  With 

iji  =  \jt  =  0 


^  +  d  ^  +/\  <3^  =  0 


(A-i) 
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Let  ua  assxsae  as  the  solixtioa  of  Sq.  (A>3)  the  following  expression: 

tdt 

OC  ^  k  a"  ;  A  =  const.  (A-t) 

Differentiating  Sq.  (A-i^)  and  substituting  into  Bq.  {A-3)  yields 

A  +  >\  +  d\+_rL  «  0  (A-5  ) 

2 

If  the  coefficients  d  andyx  were  constants,  then  two  roots  which  are 
particular  solutions  of  Eq.  {A-5)  are 

X"  -  |a|  (A-6) 

However,  since  the  coefficients  are  tiiae  variables ^  ^  constant 
does  not  satisfy  Sq.  (A-5). 

At  this  point,  let  us  nake  the  following  transfoxsation: 

y,-  a  +  1  CO  (A-T  ) 

Siibstitutlng  Sq.  (A-7)  into  Bq.  (A-5)  and  sinqtlifylng  yields 


0^  +  2J  o6)-CO^+o  +  J  CO  +  do  +  J  d.C*>+j^ 
Separating  the  real  and  imaginary  parts  of  Bq.  (A-8) 


0  a 


4  (J 
2  2CJ 


and 


{A-8) 


(A-9) 


ol 


d 


f 


(J^  -f  d  a  +  (7  (A-10  ) 

Frees  Eq.  (A-9)  it  can  be  seen  that  even  if  there  were  no  aerodynamic  dan5>iag 
of  the  re-entry  body,  i.e.,  d  equal  to  zero,  the  variable  coefficients  affect 
the  asgtlitude  of  the  oscillations  through  the  change  and  the  rate  of  change 


of  the  4aKp»d  nutural  frequency, CJ .  Coesblning  *q.  (A-9)  aal  (A-IO)  yioicU 
for 


,  ,2  -2  d  3  1  y 

o  ^  -  C  -  5 .  r  (jj)  -  5  “ 

since  the  real  part  of ^  a*  given  hy  Sq.  (A-9 )  vill  integrate 
directly,  the  eolutlon  of  Eq.  (A-3)  ie 

^  -  1/2^  dtf  jfcidt  -jlcJdt 

^  ^  =1  •  J  ♦  =2  •  J 

The  daaped  natural  frequenay,CJ>  reaaine  to  be  determined.  Before  any 
eolutlon  f or  CJ  nay  be  obtained  froo  Bq.  (A-ll),  It  is  necessary  to  kncvj^ 


(A-il) 


(A-i2) 


and  d  as  functions  of  tias. 

I^OD  an  exasdnatloQ  of  the  relative  ■agnitudes  of  the  terms  of  Bq.  (A-U), 
it  is  found  that,  at  least  for  conical  re-entry  bodies 


2  d^  dj 


(A-13) 


That  is,  the  rate  dasiylng  of  such  a  Teblcle  is  so  sjbaU  that  the  dasiped 

and  undasqped  natural  frequencies  are  for  all  practical  purposes  Identical. 

2 

fujrtboriBoro ,  from  the  definition  of_r\_  on  page  31  it  is  evident  that 


2  2 

Thus  as  a  first  approx ijsat ion  to  the  solution  of  Bq.  (A-ll) 


(A-IU) 


03  -40^  (A-l'>) 

Xquation  (A-1^)  represents  the  exact  solution  of  Bq.  (A-ll)  (still  neglecting 

2  * 

the  extreasly  small  d  and  d  terns )  if 


2  3 


-f'ST) 


D,  1 

r)  *5y- 

a  n 


(A-16) 


See  the  discussion  of  the  dcjoping  terms  in  Ref.  2. 


Also,  see  Appendix  C. 
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where  the  stibscrlpt  *1'  indicates  the  first  approxlsaation  of  the  coefficient 
2 

-A_.  If  for  the  particular  conditions  vmder  conBidaration  the  derivative 

2 

teras  of  £q.  (A-l6)  ere  not  negligible  caiiq>ared  to  CO  ,  than  a  second 

n 

iteration  snurt  be  utilised.  Thus 


f  >  I  As  3  1 

n  (o")  -  2  5- 

'  n  n 


and 


— 

-  n 

(^)'  -  («)' 
61  '  ^61^ 
‘"n  _ 

1 

di  ; ; 

n  CO 

*  2 

(A-17) 


o  2  ,,2  3 

A.  =0.  +  r 


2  2 

The  iteration  procesa  aiay  be  continued  until-T^  is  aa  close  to^  as  is 
desired.  For  the  problem  at  harkd^  the  foUoving  eoq^ression  has  been  found 
to  be  qpite  adequate: 


(A.18) 


O  .  Yh*  +  jg  k®  sin®  |« 

2 

The  tern  imder  the  radical  of  Sq.  (A-I9)  is  amch  larger  (on  the  order  of 
a  factor  of  100  at  250,000  ft  altitude)  than  the  reaaining  temi  except  at 
the  extr^ae  altitudes  at  which  the  initial  conditions  of  the  probleia  must  be 
established.  Equation  (A~12)  thus  becones 

d  dt 


(A-19) 


^ 

o 


jL  1/k  - - - 


-li 


q.*  +  ^  •in®  t 


=571  ““ 


(A-20) 


Z  -  +1  C*)dt 


wheire 
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In  obtaining  Bq.  (A-20),  it  haa  been  aasvoned  that  re-entry  starts  at  an 
aititiide  of  approxiiaately  500,000  ft.  At  this  altitude  there  is  en  angle, 
between  the  velocity  vector  and  the  longitudinal  axis  of  the  vehicle. 

The  rate  of  change  of  is  initially  zero. 

Ihe  planar  case  is  considered  in  greater  detail  in  yi3K?®^^cLix  D. 

Let  us  rettim  to  Bq.  (A-l)  and  (A-2)  and  consider  the  coupled  case  in 
^ch  the  re-entry  body  is  spinning  about  its  longitudinal  axis.  If  the 
coefficients  of  Bq.  (A-l)  and  (A-2)  were  constazrt,  solutions  might  be 
obtaiiied  by  any  one  of  several  nethods.  Perhaps  the  most  convenient  approach 
is  to  assume  an  exponential  solution  for  both  and  p  and  then  derive  the 
secular  equation  (e.g.,  see  Ref.  t).  The  secular  equation  may  be 
solved  for  the  required  roots.  An  analogous  procedure  may  be  used  for  this 
particular  problem,  even  though  the  coefficients  are  functions  of  time.* 

Let  us  assume  solutions  totxcL  and  3  of  the  following  form: 


^  =  A  e 

(A-21) 

fAdt 

3  *  B  e 

(A-22) 

where  A  and  B  are  arbitrary  constants. 

Substituting  and  p  and  their  derivatives,  as  obtained  frcci  Bq.  (A-21 ) 
and  (A-22),  into  Bq.  (A-l)  and  (A-2)  yields 

A  ^  +  d\+j^)  +  B  (X+  c)  =  0  (A-23) 

-A  R^  (X+  c )  +  B  (A  4  dA*_rL)  0  (A-24) 

#  ” 

IMs  is  true  since  the  problem  represents  the  so-called  'degenerate 
case.'  That  is,  with  zero  coupling  tbe^sC  and  3  equations  would  have 
identical  solutions. 
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vith  A  ax2d  B  Icnovn  constants,  a  non- trivial  soltztion  exists  only  if  the 
deteradnant  of  their  coefficients  equals  nsro,  thus 

2 

+  +  +K^^(X+c)^  =  o 

Ag&in  taking,^^to  be  a  coE^>lex  qjoantity 


X  =  o  +  j6) 

Substituting  Bq.  (A-26)  into  Sq.  (A-2^),  and  separating  the  real  and 


iaaginary  parts  yields 


0 


d 


CO 


26)+  R  20+  R 


20+  R 


and 


o  «  R  e  d  ,2  ,  V 

- i-y. 

^  20+  R^  20-^  R^ 


-  2 


dCO 


2  2 


20+  R  (20+  R  y 

X  X 


d  O 


20+  R 


4.  ^  : 

'  ^  » 


2RcdO  URcCa) 
-  A. 


(20+  R  f  (20^  R  (20;  R 

X  '  X  X 


R  ^ 


R  c 

+ 


CO 


(20+  R  r  ■  26>+  R  (20+  R  F  20+  R 

X  X  x'  X 

If  we  examine  the  magnitude  of  the  various  terms  of  Sq.  (A-28),  it  is 

2 

apparent  that  those  containing  c  and  d  are  very  small  eonqvared  to_A.. 
As  a  first  approximation  of  Owe  find 


(A-25) 


(A-26) 


(A-27) 


(A-28) 


(A-29) 
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An  exaaination  of  the  derivative  terns  of  Eq.  {A-28)  reveals  that  Zq. 
(A-29)  represents  a  satisfactory  apprond nation  over  an  altitude  range  of 
500,000  ft  to  sea  level.  Hitis  the  four  frequencies  of  oscillation  are: 


(A-30) 


(A-5i) 
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BqimtioQ  (A«3l},  (a)  throxigh  (d),  reveals  that  the  tvo  sets  of  roots  are  not 
ideoticalXy  daagped.  Hovever,  an  investigation  of  relative  aagnl  tudes  indicates 
that  the  middle  tvo  terms  of  Sq.  (A-3i),  (a)  and  (b),  say  be  neglected  in 
eosparlsoa  to  the  rwaalnlng  terms  in  the  range  of  altitudes  froa  2^,000  ft 
to  sea  level.  Ihtis  the  soluticms  of  Eq.  (A-l)  and  (A-2)  have  the  folloving 
form: 


cos  +  A^  jin  Z,  +  A_  cos  Z,  +  A,  sin  Z 


1  ^  ''3  3  ^  % 


(A-32) 


-  « 


0 


iA‘ 


cos  Z^  +  Bg  sin  Z^  +  coa  Z^  +  B^^  sin  Z, 


(A-33 ) 


vfaera 

Z  =fo  dt 

X  ^  i 

Z3  dt 

Of  the  eight  constants  of  integration  that  appear  in  Eq.  (A-32)  and 
(A-33),  <s3ly  four  are  independent,  since  only  foxir  total  initial  conditions 
may  be  specified.  Ihe  choice  of  Initial  cozidltions  is  dependent  upon  the 
assuaed  performance  of  the  attitxule  control  system  xdiich  positions  the  vehicle 
for  re-entry.  Thus  if  at  the  end  of  the  control  period,  the  body  has  both 
an  error  in  attitude,  and  a  residual  body  angular  rate,  then  the  follovlng 
set  of  Initial  conditions  are  reasonable: 
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a.  a, 

O  w 


a  =  a 

o  o 


=0 

If  such  a  set  of  initial  conditions  is  selected,  the  foUcnrlng  constants 
of  integration  ssay  be  obtained  frca  E^.  (A-32)  and  (A-53)  vban  the  initial 
altitude  is  such  that  the  influence  of  the  aerodynaxle  terns  is  negligible. 


A^  =  0 

(a) 

X 

(b) 

a  R 

A  ~  ®  ^  . 

(c) 

i  ■  /T 

v2 

Qn 

A  ^  V/ 

X 

(d) 

B 

1  R 

X 

(a) 

=  0 

(b) 

X 

(c) 

a  R  1/2 

g  ~  ox 

**  y: 

(d) 

(A-J-^) 


(A-3'  ) 


RM-1863-1 

ko 


For  thft  «ss\atd  eondltloas,  the  soltxtlons  of  Sq.  (A>52)  and  (A-35)  my 
he  vrltten  as  fotUotrs: 


d  dt 


a  =. 


-  H  2 

.  2  X 

U)  +  "T 

n  4 


lA 


R 


Sin 


[vV“)] 


a  R 
o  X 


1/2 


^/2 


coijz,-Z^  (o) 


aM 


■II 


t  r 

d  dt 
o 


ij 


2  V 


J  4. 

n  4 


1/4 


r 


R 

2 


jws  jz^-Zj(o)j  -  cos  Z^-Z^(o) 


a  R 
o  X 


1/2 


sin 


[^-=^3  '»>] 


(A-57) 


(A-58) 


In  ohtaialcg  the  approxlmte  solutlozis  for  Cl  and  as  given  by  Eq.  (A-37 ) 

aiBl  (A»3d)j  it  has  been  assuaaed  that  the  tena  containing  the  static  mrgin, 

2 

0  ,  is  predoadLnant  frosa  approxiaately  250,000  ft  to  sea  level.  For  the 
n 

special  ceue  la  'vdiich  the  static  laargin  is  zero,  the  four  roots  are: 


^1  »  -  Kjj  +  J  R^  (^) 

X2  •■X3  -  -  c  (A-39) 


(=) 

With  the  initial  cwaditions  given  by  Bq.  (A-34)  the  six  constants  of  integration 
are: 


RM-i8^1 
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(a) 

A 

2  R 

(b) 

(A-kO) 

X 

(c) 

aj2d 

ob- 

B  =  -  -2 

1  R 

(a) 

X 

-  0 

(b) 

{A-kl) 

B 

®3  ~  R 

(c) 

X 

For  the  case  of  zero  static  margin,  Eq. 

(A-32)  and  (A-33)  have 

solutions 

of  the  following  form: 

-  j"  c  dt 

=  e  ~  sin  R  t 

‘'x 

+  e 

0 

(A-kS) 

and 

-  f  dt  . 

P  ^  ^  ^0 

-|cat 

(A-4j) 

f  -  -  e  =~  cos  R 

"x  * 

t  4  ^  e 

X 
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APPENDIX  B 


The  Forced  Solutioaa  of  the  gomt-icms  cf  Itotlcm 


Pron  Sectioa  II,  the  coupled  equ&tlaia  inOC  and  B  are: 

•  ••p  •  ••  •  A—  (B  ) 

^  *  c  p)  =  0,,  ij^)  — ^ 


cos  t 


(B-1) 


3  *  d  P  3  -  W.  c«)  =  -  .of  (j^)  >11.  *  ( 

In  Appendix  A,  tha  hotaogeo^oxis  solxitlcsis  of  Eq,.  (B~1)  and  (B-2)  are  extunined. 
Once  the  homogeneous  solution  of  a  linear  differential  equatlco  has  been 
detera&ined,  it  is  alvays  possible  to  find  the  ccnplete  solution,  at  least  in 
on  integral  form.  In  the  present  case,  the  integrals  Involved  in  determining 
the  forced,  or  particular,  solutions  are  very  cao5>lex,  and  ae  a  result, 
appraxlmate  methods  snirt  be  utilized. 

One  of  the  first  approxiioations  that  might  be  made  is  to  neglect  the 
derivative  terms  that  occur  in  Sq.  (B-l)  and  (B-2).  An  examination  of  the 
various  forcing  functions  of  Eq.  (B-l)  and  (B-2)  indicates  that  the  forced 
or  'steady  state'  solxrtions  should  contain  primarily  very  low  frequency 
cosqxments .  Thvis 


(B-2) 


"x  =  K. 


cos  if 


{B-5) 


-  R  c  (X 


2  A-  (B+Ca) 


(B-Ji) 


SolviDg  Eq.  (B-3)  and  (B-K)  forCC  and  B  ,  and  simplifying  ylelda, 

S8  88 

respectively 


2 

-TL 


o  A_  (5  +6^  ) 


*  R  c  0 


S8 


88 


-Tl!  +  R  ^ 

X 


R 

X 


c  CO 
n 


(6  +eA  ) 

£  M 


'V*B  ^c 


sin  y 


J^L  +  R  ^ 

X 


2 

c 


(3--) 


and 


2r  :  2 


+  Rc0  +RcK^0 

X  '^SS  X  « 


68 


88 


o  2  2 

_n_  +  R  c 

X 


P  (6  +  6A  ) 

2  /  D  \  2  X 


cos  \jf 


D  2  2 

-TL  +  R  c 


(B-c) 


For  the  usxml  case  in  vhich  the  term  representing  the  static  moaaent  is 
predominant,  Eq.  (B-5)  and  (B-'j)  reduce  to  the  foUcving  form: 


CC 


K.  SiK,  ,(V6\) 


8fl  CJ  CO  2 
n  n 


cos  t 


•  (B-7) 
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and 


& 


88 


^  *  'oc) - Z - ♦ 


a 


(B-8) 


By  differaatiatlng  Eq.  (B-7)  and  (B-8)  azkL  Introducing  the  reiults  Into 
Eq.  (B-l)  and  (B-2),  a  second  approximation  for  0^  and  3  nay  be  obtained. 

88  8C 

However,  unless  the  static  sargin  of  the  vehicle  ia  very  close  to  lero, 

Bq.  (B-7)  and  (B-8)  are  excellent  representations  of  the  forced  solutions  in 

the  altitude  range  from  250,000  ft  to  sea  level. 

?roa  the  definition  of  CJ  ,  given  on  page  8  ,  and  of  given  on  page 

31  ve  find  for  and  3 
88  88 


88 


d  K 

-  dt  (y  Si  8  ®  >  (V 

»  — — — — — -  e - - - s - +  [r—rr')  - cos  f 

,  W  A 


(B-9) 


and 


B 


I  CO  g  COB  6  e 
X  X  '^ss 


-kh 


88 


Ajj  ^  (&j«-eA^) 


sin  ijr 


(B-10) 


From  Bq.  (B-9)  end  (B-IO),  it  can  be  seen  that  as  approachee  zero, 

the  steady-state  values  of  OC  and  3  approach  infinity.  But  as  A  goes  to  zero, 
2 

beccmes  very  small  and  the  original  appraximation  is  no  longer  valid 
since  the  derivative  terms  become  eigaif leant. 

For  the  special  case  in  which  A^  is  zero,  the  exact  forced  solutions 
for ^  and  3  ®ay  he  deduced  from  Eq.  (B-l)  and  (B-2).  Thus 


-J'c  dt^  Jc  dt 


Oi  «  e 

SB 


T. 


d  ^ 


ss 


(B-11) 


1^6 

and 

3  ts  0 

ss 

Upon  considering  Kq.  (B-10)  and  (B-12),  It  appears  that  as  A  becoaes  small, 

X 

Initially  the  forced  solutions  increase  in  magnitude  until  the  point  is 
reached  where  the  derivative  terms  become  important.  Prom  that  point  until 
is  zero,  the  magnitude  of  the  forced  terms  decreaae. 


(B-12) 
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APPEHDIX  C 


The  Aerodynamic  Parsunsters 


Oa  page  7  ‘tbs  lift  and  drag  forces  are  defined  In  terms  of  tvo 
parameters,  and  These  parameters,  as  indicated  by  Xq.  (sl)  and  (22), 
are  functions  of  the  sea  level  air  density,  a  body  reference  area,  and  the 
lift  curve  slope  aM  drag  coefficient,  respectively. 

Thus 


e  C>C 


(a) 


(b) 

(c) 


(Col) 


where 


\  “  2  ^o  ^ref. 

S  “  2  ^o  ^ref .  S 


The  daaqjing  coefficient,  may  be  defined  in  a  slallar  manner. 


ihus 


kh 


(C-2) 


where 


(c  +  c  ) 

a  n  •, 
<L  X 


2  I 


^o  ^ref. 


and  is  the  length  of  the  body. 


If  a  coalcal  re-entry  body  is  considered,  then  at  hypersonic  speeds 


impact  theory  yields  the  foLLoving  expressions  for  C.  ,  CL,  C  and  C 

Lt  1/  m  ] 


2  (cos^  B  -  sin^  B) 


iC 

(a) 


B  2  sin  B 


(b)  (C-3) 


m 

<1 

q- 


2 

cos  B 


"  -|(K) 

■^ref 


2  cos^  B  ■) 
'^ref 


(c) 


C  =0 
m 

dc 

O 


(d) 


where  9.^^  is  the  distance  from  the  nose  of  the  cone  to  the  center  of  gravity 
and  B  is  the  cone  half-engle.  In  terms  of  the  static  margin,  A  ,  ve  find 

n 

for  : 

'^ref 


‘^cg  2  1 

^  COB^  B 


(C-ii) 


Utilizing  Sq.  (C-3)  and  (C-4)  ve  find  for  A^,  and 


p  A  -  (cos^  D  -  sin^  B) 
o  ref 


p  A  „  sin  B 
o  ref 


0  PT  ^ 
o  ref 


0  A  A  ^  ^ 

o  X  ref  2  „ 

-  cos  B 


I8TTI  sin  B 


See  Ref.  2. 


(C-^) 


(C-t  ) 


(C-7) 
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The  loading  peLraiaater^  It,  Is  dofinad  as  follovs 


R  »  ~2S_ 
^  A 


(C-8) 


ref  D  2  A  -  sin**  B 

By  substituting  Eq.  (C-8)  into  Bq.  (C-5)“(C-7)>  reference  area  nay  be 
eliminated.  We  are  nov  in  a  position  to  evaluate  the  relative  magnitudes  of 
the  coefficients  of  Eq.  (A-l)  &7jd  (A-2).  Ihus: 
l)  The  djBnq)lng  term 


A^  kh  khri  g  2  2 

do-^Ye  +K^»Ye  1-^  (cot  B-l)  +  ■  +  -2— ^ cot  B 


72TTI  R^ 


sin  B  2  N  I 


2)  The  tndaacped  oatural  frequency 


O  O  A 


(C-9) 


^-r  “8  P  P  ^  8  .  P 

^  cot^  B  e  +  (cot^  B-l) 


.  kh  .  kh 
V  e  +  kh  V  e 


(C-10) 


P_  g  p 
^  (cot^  B-l) 


pQ  (»g)' 


TTTT 


^2 

A,  ag  p 

+  -2—2^ - cot^  B 


L72Tri  W  sia  B  2  R  I 
2 


v2 

V  e 


2kh 


The  expression  for may  be  further  simplified  by  recognising  the 
foXIorvlng  relationships: 

A^  P  hh 

Vb - idy^e  +g  sin  6 

and 


h  -  -  7  sin 


ss 
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3)  "nie  aorodynanlc  coupling  tenn 


m  ^  "  2  N 


V  e  =  (cot'"  B-l)  V  e 


(c-n) 


An  inspection  of  Eq.  (C-9)-(C-ll)  Indicates  that  for  a  cone  balf-anele 

2 

of  45  deg,  the  expressions  for  d,-rL  and  C  are  greatly  sijapiified.  IMs  is 
due  to  the  fact  that  ii^^wit  theory  yields  a  lift  curve  slc^e  equal  to  aiero 
for  such  a  case. 

2 

The  coc5>aratlve  magnitudes  of  the  three  terms  vhich  constitute may 
be  seen  in  the  following  table  for  a  typical  re-entry  cone: 


Table  C-1 


- -  ^  2H  2  2N  d  A  „  .kh 

A  .  .  0.025  n  - 

Po  «  ^  e  P  g  V  e 


Ul 


B  K  30  deg 
N  «  140 


m  B  100  slugs 

h  »  500,000  ft  4,^  X  10* 

250,000  4.69  X  10* 

50,000  4.c^  X  10* 


2.46  X  10’ 
2.46  X  lO' 


1.52  X  10 


5.57  X  10 

1.57  X  10* 


5.73  X  10 


From  Table  C-1  it  can  be  seen  that  even  with  a  static  margin  of  0.5 
2 

of  an  inch,  the  term  is,  under  the  least  favorable  circumstances,  larger 
than  the  remaining  terms  by  a  factor  on  the  order  of  100.  Ifader  the  same 
conditions,  the  dan^jing  term,  d,  and  the  coupling  term,  c,  are  roughly  one- 
tenth  the  magnitude  of  .  Thus  terms  containing  the  square  of  c  and  d, 
or  their  product,  ^Ich  occur  in  Eq.  (A-28)  may  be  neglected  as  compared  to 

.  A  negative  increase  in  sein/es  to  increase  the  disparity  in  nsagnltvideB 
between 6^^  and  the  other  terms. 
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APPffiCm  D 

PXaiiST  PvS^^gtry  AbsItsIs 

Ibe  siJziljLrlty  betve^n  the  pl«n&r  solution  obtained  In  Appendix  A  and 

« 

the  results  obtained  by  other  Investigators  is  not  isaedlately  obvious. 

By  solving  the  problen  vlth  altitude  rather  than  tlae  as  the  Independ^t 
variable,  the  degree  of  equivalence  of  Xq.  (A-20)  to  the  results  of  Bef.  2 
beccsaes  more  apparent. 

Hotlng  that 


§  =  -  V  sin  ^ 
dt  ^ss 

Xq.  (a-3  )  aay  be  traasforoed  to  the  follofviirg  form 


d^cC 

d(k 

»  AIS 

*D* 

dh*" 

dh 

a  sin  0 

«/ 

! 

1 

1 

^  sln^  0 
'^ss 

1^ 

-  '  "j  ■  '  < 


kh 


ss 


SB 


e  ■¥ 


V®  Bin®  i  “ 

'^ss 


B  sln^  ^ 


e 


2ch 


ss 


Xm  0 


By  the  method  utilized  In  Appendix  A,  ve  find 


O  (h)  S  . 


<V^)> 


sin  ^ 


kh  3 
e  ^  ^ 


cS(h: 

JtF 


ss 


2Uth) 


See  page  8,  Ref.  2. 


(1^1) 


(D-2) 


(1>«3) 
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vher«  the  second^order  terms  have  been  neglected.  As  a  first  npproxination, 


sin  0 
'^ss 

Substituting  C)j^(h)  andCJ^(h)  Into  Eq.  (D-3)  Indicates  that  Eq.  (D-*i»)  is  a 
good  approxlstation  as  long  as 


(D-4) 


_x 

I 


ss 


For  a  of  -0.5  ^  or  larger  (in  a  negative  sense),  the  inequality  of 
Eq.  (D-5)  is  valid  for  altitudes  from  250,000  ft  to  sea  level. 


If  the 


JL 


ain^  0 


term  of  Eq.  (l>-2)  is  neglected  ve  find  for  o(h): 


0(h) 


A  ^ 
*1)« 
m  sin  0 


kh 


ss 


sin  0 


88 


- 

25(hT 


Thus  for  X 


1 

2 


a 


ffl 


k  sin  0 


88 


kh 


X 


\<V*D>  .to 

1/k 

J^(h)  dh  -J^CJ(h) 


dh 


Cl  e 


+  9 


(D-5) 


(D-6) 


(D-7) 


where  0^^  is  assumed  to  be  a  constant.  A  coB5>ari8on  of  Eq.  (D-7  )  with  Sq. 
(I5e)  of  Ref.  2  reveals  that  the  terms  which  affect  the  aoqaitude  of  the 
oscillation  are  identical  except  for  a  constant  which  may  be  absorbed  in 


=1  's- 
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la  Ref.  2,  the  poiat  is  made  that  for  a  hlgh-drag  body,  it  vo'old  be 

possible  for  the  exponent  of  Eq.  (D-7)  to  be  positive,  and  thus  for  divergence 

to  occin:.  However,  an  examination  of  the  origin  of  the  drag  term  which  occurs 

in  the  eaqponent  of  the  e.ponential  in  conjunction  with  the  ijiplicit  assumptions 

that  have  been  made  by  neglecting  — s - - term  of  Eq.  (D-2),  indicates 

sin  0 

that  divergence  due  to  a  high -drag  configuration  is  not  a  serious  problem. 

As  an  intermediate  step  in  the  transformation  of  Eq.  (A-3)  to  Eq.  (D-2) 
by  the  use  of  Eq.  (D-l),  we  find,  for  the  coeJfficient  of  the  following 
expression: 


d0 


/i  dv  ^  '  ss  d  \ 

T  dh  ^  ^88  dh  "  V  sin  0  ' 

^ss 

Under  the  assumptions  of  Ref.  2,  0  is  a  constaait,  thus  /dh  is  zero,  and 

83  So 

the  speed  as  a  function  of  altitude  is  approximated  by 


^  _kh 

m  k  s  in  0  ® 


V  =  7  e  ®®  (D-8) 

o  . 

where  the  component  of  g  along  the  flight  path  has  been  neglected.  With 

dec 

these  approximations,  the  coefficient  of  ^  reduces  to  the  form  used  in 
Eq.  (D-6).  However,  the  method  of  solution  adopted  in  this  report  does  not 
require  the  above  approximations. 

IhuB 


a(h) 


1 

2 


1  ^ 
V  dh 


cot  0 


ss 


ss  dh 


\  ,  , 

a  *  \lj  ih 

sin  0 
'^88 


1  CJ(h 

2G5(h 


(D-9) 


5^^ 


and 


2  k  sin  0 


V  sin  0 
o  c 


P-/2 


jjo(li)  dh  -jfcHh)  dh’ 


+  0^6 


(D-10) 


Unless  the  lift  curve  slope  is  negative,  the  exponent  of  the  exponential 

/V  y/a 

tern  of  Eq.  (D-IO)  is  never  positive.  Tne  ratio  of  Eq.  (D-IO) 

corresponds  to  the  exponential  drag  term  of  Eq.  (D-7).  Althovigh  V  can  become 
quite  smal.l  compared  to  V^,  it  never  approaches  zero.  The  terminal  speed  of 
almost  any  conceivable  re-entry  body  vould  be  on  the  order  of  500  ft/sec. 

In  Eq.  (D-7),  the  coQqwnent  of  g  along  the  flight  path  has  been  neglected, 
and  as  a  consequence,  V  may  approach  zero,  and  the  amplitude  otcC  may  thus 
approach  infinity.  Considering  Just  the  envelope  of  the  oscillation,  as 
given  by  Eq.  (D-IO),  we  liave 


Equation  (A-20)  of  Appendix  A  is  equivalent  to  Eq.  (D-H)  if  the  initial 

p 

conditions  are  established  at  an  altitude  of  2^)0,000  ft  where  the  k  term  is 

r  2 

small,  compared  to  CJ  "  .  At  hi^r.er  altitudes,  where  CJ  becomes  very  btorI i 

n  '  '  n 


we  find  for  the  envelope: 


50  100  150  200  250  300  350  400  450 

Altitude, h  (thousands  of  ft) 

Fig.  10 — Planar  angle  of  attack  envelope  as  a  function  of  altitude 


Normalized  planer  angle  of  attack  envelope 

O  p  o 

O  ro  O)  00  o 


9^ 
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